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ABSTRACT 

We present the analysis of two long, quasi-uninterrupted RXTE observations of Cygnus X-1 that span several days within a 10 d 
interval. The spectral characteristics during this observation cover the region where previous observations have shown the source to 
be most dynamic. Despite that the source behavior on time scales of hours and days is remarkably similar to that on year time scales. 
This includes a variety of spectral/temporal correlations that previously had only been observed over CygX-l's long-term evolution. 
Furthermore, we observe a full transition from a hard to a soft spectral state that occurs within less than 2.5 hours - shorter than 
previously reported for any other similar Cyg X- 1 transition. 

We describe the spectra with a phenomenological model dominated by a broken power law, and we fit the X-ray variability power 
spectra with a combination of a cutoff power law and Lorentzian components. The spectral and timing properties are correlated: 
the power spectrum Lorentzian components have an energy-dependent amplitude, and their peak frequencies increase with photon 
spectral index. Averaged over 3.2-10 Hz, the time lag between the variability in the 4.5-5.7 keV and 9.5-15 keV bands increases with 
decreasing hardness when the variability is dominated by the Lorentzian components during the hard state. The lag is small when 
there is a large power law noise contribution, shortly after the transition to the soft state. Interestingly, the soft state not only shows 
the shortest lags, but also the longest lags when the spectrum is at its softest and faintest. We discuss our results in terms of emission 
models for black hole binaries. 

Key words. X-rays: binaries, X-rays: individual: Cyg X-1 
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2 1 ■ Introduction IWilms et al.l2006h . In the hard state, BH XRBs are generally ob- 
servable in the radi o with a flat/inverted rad io through IR spec- 

The observed behavior of most black hole (BH) X-ray bina- tmm (Fender et al . J200I l2009t iDunn et aflEoToh . X-ray vari- 

nes (XRBs) can be classified into two major states, the so- abmty aboye Ql Hz during this state . g high (rms > 2Q%) In 

called hard, power-law dominated, and the soft, thermally dom- the soft state the x . ray spectmm is dominate d by emission from 

mated state (e.g.,| Remillard & McChntock|2006i and references an accretion disk and the core radio/IR emission is quenched. 

H , therein; these authors also introduce the most extreme steep C ygX-l, however, retains a strong power law component even 

power law state seen in only a few BHs). In the hard state the • - t so ft es t X-rav soectra 
X-ray continuum can be described by a hard power law, which is 

interpreted as resulting from Comptonization in a hot (kT e » 50- In most Galactic BHs, transitions between states are com- 

lOO keV) electron plasma surrounding the accre t ion disk, or in a tnon. The state behavior generally can be described as part of 

iet (iDove et alJ[l99 8: Haardt & Maraschil [19931: IZdziarski et alJ a characteristic evolution of spectral and timing states that de- 
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2002; Marko ffet al J 120051) . A strong contribu tion of jet syn- pends on the previous source history. Over the outburst of an 

chrotron emission is found in some sour ces (Markoff et al. 2001; X-ray transient, the source is seen to move along a q-shaped 

iRussell et aT1l2010t lLaurentetafll201lh . In the case of CygX-l track in a hardness intensity d iagram (HID; e.g.. iFender et alJ 

the X-ray continuum can be modeled with a broken power 2004; Homan & Belloni 2005). The outburst starts in the hard 

law hardening above a break energy E b[e . dk a lOkeV and a state; during the initial phase, the source brightens without sig- 

high-energy cutoff with a folding energy of ~100keV (e.g., nificantly changing its spectral shape, and the radio/IR flux in- 

creases. At some point, the X-ray spectrum starts to soften. 

Send offprint requests to: M. Bock, The radio/IR flux becomes more erratic. In some s ources radio 

e-mail: Moritz.Boeck@sternwarte.uni-erlangen.de ejections were seen during this intermediate state (Corbel et al. 
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12004 iFenderet al.1 [2004 120061: IWilms et alJl2007h . After the 
state transition, the source spectrum is soft and no core radio/IR 
emission is detected (Fender et al. 2009). During this soft state 
phase the bolometric source flux decreases, until at some lower 
luminosity the X-ray spectrum hardens, and the source becomes 
detectable in the radio/IR again. During the final part of the 
outburst the source is back in the hard state and the source 
luminosity decreases until th e XRB ceases to be observable 
dFender et a l. 2004; Belloni et al. 2006). 

The observed behavior of BH XRBs is determined by the in- 
terplay between emission from the accretion disk and from the 
regions of the accretion inflow/outflow in which the hard X-rays 
and the radio emission are generated. Transitions between the 
two states are therefore especially suited to place constraints on 
the physics of these components. Despite the fact that the out- 
bursts can last weeks or months, the transitions bet ween hard and 
soft st ates are for most sources rare and often fast. Belloni et al. 
(2006), for example, report a hard to soft transition of GX 339-4 
that occurred in less than 10 h. For this reason, the behavior dur- 
ing state transitions is not well studied. In addition, comparisons 
between different sources are complicated because that the influ- 
ence of other variables, such as the variation of the mass accre- 
tion rate, is difficult to ascertain. Here we will present a detailed 
analysis of a state transition as observed in a single object. 

The canonical BH candidate Cygnus X-l is a persistent XRB 
that spends most of the time in the h ard state; however, it shows 
frequent excursions into soft states dWilms et alj|2006l and ref- 
erences t herein), including aborted trans itions ("failed state tran- 
sitions"; iPottschmidt et alj|2000l 120031) . Because of its bright- 
ness and persistent emission, CygX-1 is one of the best stud- 
ied BHs. In addition to studies o f snapshot obs ervations (e.g. , 
Gierlihsk i et alj[l999 : Di Sal vo et al Il200it iFrontera et alfc OQl; 
Makishim a et alj20 08: Hanke et al. 2009), most recent work has 
concentrated on the long-term behavior of the source from a 
monitoring campaign with the Rossi X-Ray Timing Explorer 
(RXTE) at a two w eek sampling tha t started in 1998 and is 
still ongoing ( Pottschmidt et al. 2003; iGleissner et alJl2004alfb1: 
lAxelsson e t al. 2005; Wilms et al. 2006, and references therein). 
The motivation of the present work was to fill the gap between 
the biweekly time scale and individual RXTE pointings with an 
intensive quasi-continuous observation spanning 10 days. Our 
aim was to determine differences and similarities between the 
behavior of Cyg X- 1 on time scales of minutes, hours, days, and 
the long-term behavior found during the monitoring campaign. 
As we will show in the following, this observation covered an 
especially interesting time when a state transition occurred. 

The remainder of this paper is structured as follows: In 
Sect. |2]we describe the observations and data analysis methods, 
the results of which are presented in Sects. [3] and [4] The results 
are discussed in Sect. [5] 

2. Observations and data analysis 

This paper is based on 78 RXTE orbits of observations of 
Cygnus X-l on 2005 February 01-11 (MJD 53402-53412; 
RXTE proposal ID 90127). We used data from all three 
instruments onboard RXTE, namely the All Sky Monitor 
(ASM; Levin eet al 1 19961) . the Proportional Counter Array 
(PCA; iJahoda et all 2006) an d the High Energy X-ray Timing 
Experiment (HEXTE~lRothschil cTetaLl [l998). We reduced 
data from PCA an d HEXTE with HEASOFT 6.3.1 following 
Wilms et al. (2006), but using data from all xenon layers and 
making use of the improvement of background models, which 
now allow us to consider data taken at least 10 minutes after pas- 



Table 1. Energy bands of the RXTE -PCA used in this paper 
according to the energy conversion of PCA Epoch 5. 

PCA Channels Energy (PCU0) Energy (PCU1234) 



11-13 

23-35 



4.53-5.82 keV 
9.72-15.4 keV 



4.49-5.71 keV 
9.40-14.76 keV 



sages through the South Atlantic Anomaly (SAA; see lFtirst et alj 
2009). Spectra and lightcurves with a time resolution of 16 s 
were obtained using the standard2f mode PCA data. To charac- 
terize the source's spectral shape on this short time scale using 
hardness studies, light curves were extracted for the two chan- 
nels listed in Table Q] as well as for the whole energy range of 
the PCA. For the timing analysis light curves with a time resolu- 
tion of ~2 ms (2~ 9 s) were extracted in the same energy channels. 
The calculation of Fourier-frequency dependent quantities such 
as pow er spectra, X-ray tim e lags, and the coherence function 
follows iNowak et all (Il999al) . Data modelin g and all other data 
reduction was performed with ISIS 1.4.9 (Houck& Denicola 
2000; Nob le et alj |2006; Nobl e & Nowak||2008l) . 

To study the spectral variability of the source, we extracted 
spectra for each RXTE orbit, yielding 78 spectra with a PCA 
exposure of about 3 ks each. We considered PCA data between 
2.8 keV and 50keV, and HEXTE data in the energy range 18- 
250keV. For the fourth and fifth PCA bin (2.8-3.2keV and 3.2- 
3.6 keV) we added systematic errors of 1% and 0.5%, respec- 
tivelyQ A multiplicative constant was used to account for the 
different PCA and HEXTE flux calibration. The constant for the 
HEXTE had typical values of about 85% with respect to the 
PCA. 

Hard X-ray spectra of BHs are traditionally explained us- 
ing Comptonization, either in a Comp ton corona or in a jet 
dWilms et al.l l2006: Mar koff et aIll2005L and references therein). 
A disadvantage of these models, however, is that the spectral 
parameters obtained from these m odels are often not unique 
dNowakl 120081 lNowaketai1l201ll) . It is still unclear which of 
the several models used in the literatu re is the correct physi- 
cal interpretation of the data. Following Wilms et al. (2006), we 
therefore modeled the X-ray spectra using a purely empirical 
spectral model, an absorbed broken power law with exponential 
cutoff. In the following, we will denote the power law indices 
below and above the break energy £break ~ 10keV with Fi and 
T2, respectively. A Gaussian profile was added to model the Fe 
Kff-line at 6.4 keV. For the thermal emissi on of the accretion 
disk the diskbb model (Mitsuda etal.ll 19841) was applied. This 
phenomenological model with Fi > F2 describes the continuum 
with its spectral hardening above 10 keV very well, with reduced 
X 2 values around 1 . We note that it is possible to map a spectrum 
defined through Fi and F2 to the mo r e phy sical models, using 
the correlations found bv l Wilms et alJ {2006). Clear correlations 
between Fi and the other parameters of the empirical model al- 
low us to identify the X-ray spectral state of the source with just 
one parameter, simplifying the analysis of the state dependence 
of X-ray timing and radio properties. 



1 Based on modeling all Cyg X- 1 spectra from the years-long RXTE 
monitoring campaign we found this selection for the two low-energy 
bins to provide reliable absorption and disk parameters. We did not add 
the recommended systematic errors of 0.5% to the remaining PCA bins 
because the reduced^ 2 of fits to hard state observations of Cyg X-l can 
be significantly below 1 in this case, indicating an overestimation of the 
uncertainties. 
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Fig. 1. Left: Hardness-intensity-diagram of CygX-1 using RXTE- ASM data. The hardness is defined here as the ratio of the 5- 
12keV and the 1.5-3 keV ASM count rates; the intensity is the total 1.5-12keV ASM count rate (see also iFender et al.l I2006L 
their Fig. 5). The grayscale represents the single-dwell (90s) ASM-data from 1996 January until 2010 August. Their characteristic 
density forms two regions in the HID, indicating two different states. The data analyzed in this work (black symbols) cover the 
transition region. Right: The transitional nature of the observation is confirmed by the HID calculated from 16 s PCA data. 



In addition to the X-ray data, quasi-simultaneous radio 
observations were made at 15 GHz with the University of 
Cambridge's Ryle Telescope. The time resolution was 8 s, using 
interleaved observations of a nearby phase reference source and 
observations of 3C48 and 3C286 as flux density calibrators. The 
telescope uses linearly-polarized feeds, and the measurements 
represent Stokes' I+Q; the fl ux density s cale i s believed to be 
consistent with the scale of Baar s et all dl977h . after allowing 
for the difference in polarizations. The overall amplitude scale 
is believed to be accurate to about 5%, and the noise level on a 
single 8-s sample is about 12mJy (the data shown in Fig.|2]have 
been averaged over lOmin, and so the thermal noise is about 
1.4mJy). 



3. X-ray spectral variability 

3.1. An hour-long state transition 

To characterize the evolution of the source's overall spectral 
properties, in Fig. Q]we show the HID of CygX-1 as obtained 
from the single-dwell 90 s i?XTE-ASM data taken between 1996 
and 2010. In this diagram, the hard state is found at a hardness of 
>1 and at count rates <40cps. Transitions into the soft state are 
characterized in the ASM by a softening of the spectrum corre- 
lated with a brightness increase, followed by a further softening 
and a brightness decrease, leading to a A-like shape in this HID. 
The right panel of Fig. Q] shows a similar HID, calculated from 
the RXTE-PCA 16 s lightcurves during the 2005 February ob- 
servation only. The figure shows that at this time resolution a 
A-shaped transitional region between hard and soft state in the 
HID is covered. In the following we will consider the two re- 
gions covered in the HID obtained with .RXTE-PCA as a "hard" 
and "soft" region instead of "hard intermediate" and "soft inter- 
mediate" regions. 



In order to investigate the potential occurrence of a state tran- 
sition in more detail we now turn to the spectral parameters of 
the broken power-law modeling of our orbit-resolved PCA and 
HEXTE spectra. The X-ray and radio lightcurves and the evolu- 
tion of the spectral parameters are shown in Fig. [2] 

During the 10 days of the 2005 February obser vation, we 
find 2 .0 < Fi < 2.7. For broken power-law models, Wilms etalj 
(2006) show that F; > 2.1 is consistent with transitional and 
soft states. Thus the range of Ti-values during our observa- 
tion confirms the conclusion from the HID that CygX-1 spent 
these days close to the transitional region between the hard 
and the soft state. A closer inspection of Fig. [2] shows that at 
MJD 53409.7, from one RXTE pointing to the next, CygX-1 
changed its photon index from Fi=2.39 to Ti=2.55 within only 
2.25 hours (Fig. [2] panel e) and moved in the HID from the 
harder region to the softer region. No other transition between 
the two regions was included in these data. The change is also 
seen in the more prominent disk, whose peak temperature in- 
creased to 0.25-0. 40 keV and whose inner radius moved inward 
to a mean value of 17 ± 3 re from 56 + 5 re before the transitionQ 
A drop in Nu after the transition could be caused by increased 
ionization of the stellar wind by the much softer X-ray spectrum. 
The duration of the transition was only a few hours, much shorter 
than the day-long transitio ns previously reported for CygX-1 
(e.g. JCui et al.lfl997l 120021) and comparable with th e very short 
transition seen in GX 339-4 bv lBelloni etail d2006l) . 



2 Values of n„ were obtained from the diskbb normalization assum- 
ing a d istance of 2 kpc, a b lack hole mass of 10 M , and an inclination 
of 35° dHerreroet all 19951) . 
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Fig. 2. Temporal evolution of the count rates and spectral parameters of CygX-1 during the observation period. The upper three 
panels are the light curves of Cyg X-l measured with (a) RXTE- ASM (1.5-12keV), (b) PCA (total energy range 2-60 keV) and (c) 
the Ryle Telescope (at 15 GHz), (d) The hardness is the same as in Fig. Q~2] namely the ratio of the PCA count rates in the energy 
bands (9.5-15 keV)/(4. 5-5. 8 keV). The bottom panels show the temporal evolution of (e) the photon index F\, (f) the photon index 
above the spectral break T2, (g) the cutoff and (h) folding energy of the high-energy cutoff, (i) the hydrogen column, (j) the inner 
temperature and (k) the inner radius of the accretion disk, and (1) the reduced^- 2 of the best fits. The dashed vertical line indicates 
the transition. 
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Photon Index I\ (E < E bleak ) 

Fig. 3. Correlation between the photon indices below and above 
the break energy, derived here for single orbit observations, is 
identical to the linear correlation present in the long -term spec- 
tral evolution from 1999-2004 (IWilms et al.ll2006l) . Note that 
this correlation extends across hard and soft states. 



3.2. Spectral variability on time scales from hours to days 

In contrast to previous observations, which consisted either of 
a few hour-long pointed observations or of several ks long ob- 
servations spaced by weeks, the present data set also allows us 
to study the spectral variability of CygX-1 during a transitional 
state. Figure [3] compares the variability of the two photon in- 
dices, Fi and F2, during the 2005 February observation, with 
that observed in the long-term monitoring between 1999 and 
2004 (Wilms et al. 2006). During the 10 d of our observation, 
the variation of Ti and F2 covered ~33% of the parameter space 
covered by the long-term monitoring, where 1.6 < Fi < 3.4 was 
found. 

This general similarity between the monitoring on time 
scales of hours and weeks also extends to the broad-band radio to 
X-ray spectral variability. The relation between the RXTE orbit- 
wise determined X-ray photon index F\ and the simultaneously 
measured 15 GHz radio flux, was consistent with the relation 
yielded by complete observations in the long-term monitoring 
campaign (Fig. |4). In the hourly and in the weekly monitoring 
we observe that the radio flux is seen to increase while the X-ray 
spectrum softens up to Fi a 2.2, and the radio lightcurve shows 
flaring behavior. Above this threshold, i.e., during the transition 
into the soft state, the radio flux decreases. 

These results indicate that the overall spectral variability of 
CygX-1 on time scales of hours is similar to that seen on time 
scales of months. The large scatter seen from one monitoring ob- 
servation to the next is due to the fact that the variability of the 
properties of the source on time scales of hours covers a signif- 
icant fraction of the overall variability of the source seen in the 
X-ray monitoring. Because the radio and X-ray fluxes also show 
similar behavior down to time scales of hours, the region where 
the radio and X-ray flux are produced must be small. 




2.5 3 

Photon Index T 1 



3.5 



Fig. 4. Relation between the 15 GHz radio flux and the photon 
index Ti . 



We expect the tight correlation between the X-rays and 
the radio to break down on time scales of minutes, es- 
pecially during flaring episodes, which are interpreted as 
ejections of hot electron bubb l es from the accretion flow 
into the jet dvanderLaanlll966t iHjellming & Johnston! 119881: 



iGleissner et al.l l2004al) . Because the radio emission is pro- 
duced farther downstream in the jet, a delay between the ra- 
dio and X-rays is expected during th ese events, a s was in- 
deed seen in the one flare reported by IWilms et al.l (120071) for 
Cyg X-l and is also observed in ot her microquasa rs, in particular 
in GRS 1915+1 05 (iMirabelet all ' 1998: Kl ein- Wolt et alJl200l: 
iPrat et al]l2010l). where it was interpreted as the ejection of the 
corona (Rodri guez et al.l 120081) . The most likely flaring event 
during our observation occurred directly after the state transi- 
tion, when the PCA light curve in its total energy range of about 
2-60 keV peaked to 5000 counts PCU , one of the highest mea- 
sured Cyg X- 1 count rates ever observed with RXTE . This count 
rate was reached although the absorption was high with a hy- 
drogen column of 3.6*2»10 22 crrT 2 . Unfortunately, however, no 
simultaneous radio data are available during this time. 



4. X-ray timing behavior 

Having discussed the spectral variability on time scales of ks, 
we now turn to the behavior of Cyg X-l on time scales of sec- 
onds. We will first characterize this behavior through standard 
power spectrum analysis in several energy bands, and then dis- 
cuss the variability through Fourier cross-spectral analysis tech- 
niques, namely X-ray time lags and the coherence function. 



4.1. Power spectra 

The calculation o f the power spectral d ensit y (PSD) was per- 
formed following iNowak et; a l. (1999aj) and IPottschmidt et al.l 
(120031) . Power spectra are normalized according to 
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o.i 1 

Frequency / [Hz 

Fig. 5. Examples of characteristic power spectra. The solid gray lines show the Lorentzian profiles, and the dashed gray lines a cutoff 
power law Left: PSD from an observation with one of the hardest spectra in our data set. /xPSD(/) vanishes for low frequencies 
and shows a double-humped shape, with no requirement for a power law in addition to the two Lorentzians. There is a slight 
indication for a third Lorentzian profile with a peak frequency between about 40-80 Hz. Middle: PSD obtained from a transitional 
state, requiring an additional power law at low frequencies. Right: PSD in the soft state, with yet stronger contribution of the power 
law. This sequence of power spectra shows that the Lorentzian peaks shift to higher frequencies as the X-ray spectrum becomes 
softer. 



lMivamotoetal.l(ll99lh . 



i.e. 



JPSD(/) d/ = I-) 



(1) 



where fi is the mean value of the light curve and <x its standard 
deviation. The PSDs were modeled as the sum of a power law 
with an exponential cutoff and two Lorentzian profiles L\ and L^ . 
We parameterize L, by the peak frequency v, where / x L,-(/) 
reaches its maximum, the quality factor Qj = /j/A/^fwhM; an d 
the strength rms,, where rms, = J L,-(/) d/ corresponds to the 
contribution to the relative root mean square variability. Using 
the peak frequency as independent parameter instead of the res- 
onance frequency 



fi 



2&v*(l+4G?)' 



-1/2 



(2) 



facilitates the fitting procedure, because v,- is much better con- 
strained and less dependent on Q, than fi, especially for broad 
Lorentzians with low Q,. v, is also the re levant parameter for cor- 
relations (e.g., Pottschm idt et alJ f2003. and references therein). 
L[ has thus the form 



Li(f) = A ! (vf-2/ i / + / 2 )" 
where /; is defined by Eq. (|2]l and 



At = rms?V; (l+4g?) 



-1/2 



- + arctan(2gj) 



(3) 



(4) 



This parameterization naturally comprises "zero frequency- 
centered Lorentzians" (e.g., Nowak 2000) for Qj - (and thus 
fi = 0), but Vi > 0. 

In our earlier work based on data from 1998 to 2001, we 
requir ed up to four Lorentzian s to model the typical power spec- 
trum (Pottschmidt et al. 2003). For the data analyzed here, two 
broad Lorentzians in addition to a power law with an expo- 
nential cutoff describ e the PSDs accurately enough (see also 
lAxelsson et al.l 120051) . The major reason for this difference is 




Peak frequency u^ [Hz] 

Fig. 6. Correlation between the peak frequencies v\ and V2 of the 
two Lorentzians. The gray line is the linear fit of Table[2] 

Table 2. Corre lation coefficients p an d weighted linear fits using 
the method of Fasano & Viol d!988l) between the Lorentzians' 
peak frequencies v, and the photon index Fj (valid for 2.0 < 
Fi < 2.7). 

p linear correlation 

0.98 ln(vi/Hz) = (-8.18 ± 0.02) + (3.70 ± 0.01) ■ T, 

0.96 ln(v 2 /Hz) = (-2.03 ± 0.01) + (1.70 ± 0.01) • T, 

0.98 ln(v,/Hz) = (1.743 ± 0.007) + (0.489 ± 0.017) • ln(v 2 /Hz) 



u 
a 

= 
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Fig. 7. Variation of the 4.5-5.8 keV power spectrum with the photon index T\. A darker color code indicates a higher value of 
/ x PSD(/). Dashes at the bottom show the photon indices of the individual spectra. Intermediate values are interpolated with a 
Gaussian blur. With increasing photon index the variability components shift to higher frequencies and their intensity decreases (in 
this energy band L2 completely fades out toward H ~ 2.6). 
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Fig. 8. Correlation between the photon index and the peak fre- 
quencies of the Lorentzian profiles. The gray lines are the linear 
fits of Table 
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Fig. 9. PSDs of the observation 90127-01-01-05 in the low 
(open squares) and high (black diamonds) energy band. The pho- 
ton index of this observation is Ti = 2. 33^'^. 



that the 1998-2001 data include observations with harder spectra 
than the data analyzed here (see also below and Fig. [5]). Because 
a strongly rising power law with strong cutoff can affect the rms 
of broad Lorentzians, we froze the index of the cutoff power 
law to 1 . In combination with the Lorentzian profiles, this 1 // 



"flicker noise" describes the data well, but it does not extend to 
the highest observed frequencies and a cutoff to the power law 
was required. 
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4.1 .1 . Correlation with the X-ray spectrum 

During the observation the peak frequency of L\, v\, varied be- 
tween 0.5 and 6 Hz, while the peak frequency of Li varied be- 
tween 4 and 30 Hz. These peak frequencies are clearly correlated 
(Fig. [6}; interestingly, we find v\ oz yfv^ (see Table |2). Figure [5] 
displays typical power spectra from the hardest to the softest ob- 
servations. The shape of the power spectra clearly depends on 
the spectral state, and exhibits the following features: 

- The PSDs of observations with a hard X-ray spectrum show 
a clear double humped profile. Individual Lorentzians can be 
clearly identified. A few of the hardest observations in this 
sample (T\ < 2.1) show slight indications of the presence of 
a third broad Lorentzian with a peak frequency in the range 
of 30-60 Hz. The contribution of the low frequencies to the 
variability of the corresponding light curve is very small, i.e., 
/ X PSD(/) decreases strongly toward low frequencies and 
no power law component is apparent. 

- With softening X-ray spectrum, the low- frequency (< 
0.1 Hz) variability increases. In this frequency range the PSD 
has a power law shape. The Lorentzian components shift to 
higher frequencies. 

- In the spectrally softest observations the power law with ex- 
ponential cutoff dominates the PSD. The rms of the second 
Lorentzian Li, which is located at higher frequencies, is sig- 
nificantly reduced with respect to L\ . 

Figure [7j shows how the power spectrum changes as a func- 
tion of the photon index Fi. The image shows the humps of 
the PSDs as dark bands, which become weaker when the X-ray 
spectrum gets softer. In addition it is apparent that the peak fre- 
quencies of the Lorentzians are clearly correlated with the pho- 
ton index T\ (see also Fig. [8}. See Table|2]for the numerical rela- 
tion between the Lorentzian parameters and the spectral shape. 



4.1 .2. Energy dependence of the variability 

Figure [9] shows the shapes of two power spectra calculated for 
the same observation but in different energy bands. While the 
qualitative shape of the PSDs is similar and the corresponding 
peak frequencies are found to agree within their uncertainties, 
the strength of the individual components of the PSD is clearly 
energy-dependent. The most obvious disparity is the strength 
of the humps. This energy dependency can be quantified using 
the Lorentzian model. The peak frequencies and quality fac- 
tors of the individual Lorentzians are generally energy-band- 
independent within their statistical uncertainty. The rms of the 
Lorentzians, however, is strongly energy-band-dependent. 

We fitted the power spectra from both energy bands simul- 
taneously, constraining v, and Qj to have the same value in all 
energy bands and leaving the Lorentzian rms a function of the 
energy band. This approach allows us to localize the Lorentzians 
in the few cases when they are almost faded out in one energy 
band but still visible in the other one. Figure[10]shows the energy 
dependence of the rms for both Lorentzians. 

The strength of L\ in the 9.5-15 keV energy band is only 
about two thirds of that in the 4.5-5.8 keV band. With an rms 
variability of >10%, L\ was visible in all the low-energy data 
analyzed here. The highest rms values for L^ are found in the 
observations with the hardest X-ray spectra (see Figs. 171 andfTUb. 
At these high values, the rms for L^ is greater in the 4.5-5.8 keV 
band than in the 9.5-15 keV band. As the spectrum softens, rms2 
decreases in both energy bands, but decreases more significantly 
in the 4.5-5.8 keV energy band. In the softest observations, the 



L2 component cannot be detected in the 4.5-5.8keV band, while 
it remains visible in the 9.5-15 keV energy band with an rms 
value of a few per cent (Fig. [10] right). In the soft state a strong 
power law component is observed at the frequencies where the 
L\ component dominates during the hard state. 

4.2. Time lags and coherence function 

While light curves at different energies are generally corre- 
lated, this relationship is frequency-dependent. The degree of 
linear correlation between the two light curves can be mea- 
sured with the coherence function, which depends on the Fourier 
freque ncies (a detailed discussion is given by iNowak et al.l 
1999a, § 3.1). For two correlated time series one can de- 
fine a Fou rier frequency-dependent delay or time lag 6,{fj) 
(Nowak & Vaughan 1996) as 



Wj) 



Kg[S*(fj)H(fj)] ®H(fj)~®S(fj) 



(5) 



2nfj 2nfj 

Here S(fj) and H(fj) are the discrete Fourier transformations 
of the soft and hard X-ray lightcurves s(t,) and h(ti), and <&$(/,) 
and <t>#(/ ) are the corresponding phases. In Eq.|5] the phase was 
chosen such that a positive time lag means that the light curve 
in the hard energy band is delayed with respect to the soft en- 
ergy band, as is common for BH XRBs. A detailed description 
of the calculation of time lags and coherence, including their in- 
terpretation and the calculation of their uncertainties, is given by 
INowak & Vaughanl d 1996b . 

The lag is generally a strong function of Fourier frequency. 
To obtain a simpler quantity, we averaged the time lag spec- 
trum and the coherenc e function in the 3.2-10 Hz band (see 
iPottschmidt et al.ll20 03). We use these averaged coherence and 
time lag values in the following. Their evolution is shown in 
Fig-QTJ In this data set the coherence was close to 1 and the time 
lag was usually longer than ~4 ms. After the transition to the soft 
state there are observations in which the time lag drops to about 
2 ms. The time lags vary strongly between about 2 and 1 1 ms af- 
ter the transition, and the coherence changes rapidly. This behav- 
ior is flux-dependent: While short time lags occur at the highest 
count rates, the time lags are long during the soft state observa- 
tions with relatively low count rates. 

Figure Q~2] again shows the HID of the present observation. 
The X-ray hardness and intensity in Fig. [12] are calculated on a 
resolution of 16 s, while the X-ray lags, indicated in the figure 
by grayscale, have been calculated on a time scale of individual 
RXTE orbits, which typically had an exposure of 3 ks. Different 
regions in the HID are dominated by different characteristic lag 
values. In the hard state the time lag is shortest at the largest 
hardness and increases as the X-ray spectrum softens. In the soft 
state the X-ray time lag is short at the highest count rates, while 
it is strongly enhanced at low count rates. 

The shortest hard st ate lag found (about 4 m s) is above the 
lowest values found by Pottschmidt et al. (2003), which is con- 
sistent with the fact that there are only observations with Fi > 
2.0 in our data set. Observations with harder X-ray spectra and 
shorter lags are not included in the sample analyzed here. 

Previous work indicated that the time lags in typical soft 
states were quite short (2-3 ms; comparable to the shortest lags 
seen in hard states), and that the time la gs increased during tran- 
sitional states (Pottschmidt et al. 2003). The temporal evolution 
shown in Fig. QT] reveals that the time lags are short (2-3 ms) 
when the variability contribution of the power law noise is large, 
whereas the Lorentzian variabilty close to the transition is as- 
sociated with long time lags (4-1 1 ms depending on spectral 
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Fig. 10. Strength of the Lorentzian profiles in the two energy bands, rms,j,igh vs. rms,i ow , for L\ (left) and L2 (right). The dashed 
line is the identity. 



hardne ss). This is consistent with the results of lPottschmidt et all 
(2003), where the shortest time lags in the soft state arose from 
variability power spectra dominated by a cutoff power law, as 
opposed to Lorentzian shape. 



5. Discussion and summary 

5.1. Discussion 

5.1.1. State transition 

The changes of the emission properties of Cyg X-l are continu- 
ous. Because the source does not follow the complete q-shaped 
track in a hardness inten sity diagram, which is found in out- 
bursts of X-ray transient s dFender et al.ll2004l : lHoman & Bellonil 
I2005t iDunnetalJbolOl) , it is difficult to specify well-defined 
borders between its different emission states. We find the follow- 
ing changes of properties: In the classical hard state of Cyg X-l 
strong radio emission is observed, which is thought to originate 
from a steady jet. The intensity of the radio emission increases 
as the spectrum softens up to about Fi » 2. 1. Above this value , 
when Cyg X-l leaves its classical hard state (Wil ms et alj |2006). 
the radio flux decrease s in correlation with the hard X-ray flux 
dGleissner et alJl2004ah . 

In the transitional state data analyzed here, we see an ad- 
ditional change of the emission properties at F; w 2.5, where 
the HID separates and the behavior of the time lags changes. 
Remarkably, this result is indicated by both spectral and tim- 
ing analysis. The data points fall into two clearly separated re- 
gions in the HID, also corresponding to different power spectral 
shapes. Transient BH XRBs only drastically change their tim- 
ing properties during transitions from "hard intermediate states" 
to "soft intermediate states". The fractional rms drops because 
of the disappearance of flat-top noise components (described as 
Lorentzians in this paper) in favor of a 1 // (power law) noise 
(Miyamoto et al., 1993; Takizawa et al., 1997; Nespoli et al., 
2003; Homan & Belloni, 2005; Casella et al., 2004; Casella 



et al., 2005; Belloni et al., 2005; Belloni, 2009). In addition, 
the "soft intermediate state" power spectra may contain specific 
QPOs referred to as type-A and type-B (Casella et al., 2005). 

iNesporietai] (120031) found a QPO at about 6 Hz in an ob- 
servation of GX 339-4 during a transition from the hard to the 
soft state. The frequency of the QPO was correlated to the PCA 
count rate. The power spectra on short time scales (16 s) of 
the data analyzed here did not show a comparable QPO, but 
it indicated that the peak frequencies of the Lorentzian pro- 
files are correlated with the PCA count rate even on time scales 
shorter than one RXTE orbit. During the hard state the PCA 
count rate is anticorrelated with the spectral hardness. Similar 
to CygX-1, GX 339-4 shows increasing time lags as the soft 
count rate increases in the hard state and the flat-top noise disap- 
pears when the GX 3 39-4 softens into a "soft intermediate state" 
ilBelloni et afl2005l) . 



5.1.2. Noise components 

Our model consisting of only two Lorentzian profiles in addi- 
tion to a power law with exponential cutoff is clearly justified 
by the structure of the power spectra shown in Fig. [7] A third 
Lorentzian is only slightly in dicated in a few of our hardest ob- 
servations (consistent with Pottschmid t et alJ l 2003h . but is not 
required in gener al. Our ident ification of the Lorentzian profiles 
is equal to that of lAxelsson et al.l (2005, see also the discussion 
in their Sect. 4.1), and our results are consistent with theirs: The 
V1-V2 relation is identical, and it is also visible here that rms2 
decreases uniformly with the peak frequency V2- 

Power spectra in the two energy bands used here can be mod- 
eled simultaneously with the normalization as the only energy- 
dependent parameter. This means that a single physical process 
is responsible for producing the individual variability compo- 
nents, but that the contribution of each process to the variability 
is energy-dependent. 
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Fig. 11. Evolution of the timing parameters during the observation. The upper panel (a) shows the PCA light curve of Fig.|2j3 for 
comparison, (b) the averaged time lag in the 3-10 Hz band, (c) coherence in the 3-10 Hz band, and (d) Lorentzian peak frequencies. 
The following panels show the rms in both energy bands (circles for the 4.5-5.8 keV band and crosses for the 9.5-15 keV band) of 
the components: (e) rms of L\ , (f) of L2, (g) the power law, and (h) the total rms. 



There are different physical models describing the origin 
and the behavior of the noise components observed in the 
power spectra. Some models connect the observed frequen- 
cies with fundamental general relativi stic ones of p erturbed or- 
bits near the compact object (see iNowak & Lehrt [19981 for a 
summary). Among others, mechanisms discussed are perturbed 



epicyclic frequencies, the Lense-Thirring precession frequency, 
and "diskoseismic" modes that are related to maximum epicyclic 
frequency. Our observed peak frequencies, 0.5-30 Hz, corre- 
spond to epicyclic frequencies at about 22-345 gravita tional 
radii for a 1OM compact object. INowak & Leh r (1998) point 
out that the "diskoseismic" frequencies are very sensitive to lu- 
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minosity changes, have a low rms variability, and thus cannot 
be applied to systems with >10% rms as in our observation. In 
addition, it is difficult to explain the spectral hardness of the vari- 
ability components. Compared to the Lorentzian profile L\, the 
one at high frequencies, Li, is more dominant in the high-energy 
band, which is clearly above energies at which the disk emis- 
sion contributes strongly, thus it is unlikely that both variabilities 
originate from the same region. Wilkinson & Uttlev (2009) ob- 
serve a low- frequency Lo r entzia n profile originating from the ac- 
cretion disk. lUttlev et all d201lh point out that the disk can pro- 



duce high-frequency variations as well, which would not show 
up as strong soft variations, provided that they occur at small 
disk radii so that only a small fraction of disk flux is modulated 
as they propagate inward. In contrast to L\ , the absolute variabil- 
ity contributions of La, and also that of the 1 // noise component, 
are comparable in both energ y bands, however. 

iPsaltis & Normanl yOOO) postulate a transition radius in the 
inner disk acting as a low band-pass filter, which yields a dynam- 
ical model for the production of modulations, with a strong res- 
onance to orbital, periastron- and nodal-precession. The model 
accounts for the observed PSDs and explains correlations be- 
tween the characteristic frequencies. We find v\ °c y/v2, which 
cannot be identified clea rly with one of the correlations between 
the modes calculated by Psaltis & Norman (2000). An extrapola- 
tion of this relation implies a characteristic frequency of ~30 Hz, 
at which vj = v%. 

In the soft state the variability (rms) at high frequencies is 
lower, indicating a relatively steady state within the analyzed 
time scales. A possible interpretation of the power law arising 
in the PSDs as the spectrum softens is flicker noise in the ac- 
cretion disk dLvubarskiill 1 9971) . The data analyzed here indicate, 
however, that the variability contribution of the 1 // component 
is strong in both energy bands and slightly stronger in the high- 
energy band compared to the low-energy band in the cases of 
strong 1// noise contribution, i.e., when the X-ray spectrum is 
soft (Fig. fT3l . A detaile d interpretation o f the s oft state power 
spectra is presented by Misra & Zdziar ski (2008), who consider 
linear and no n- linear responses of a damped harmonic oscillator, 
which is assumed to govern the disk variability at each radius. 

The dependence of the Lorentzians' rms on their peak fre- 
quency (Fig. [15]) is similar to that presented by Pottschmid t et al.l 
(2003) in Figure 6 of their work. There is, however, one slight 
difference: Their Figure 6 indicates a unique relation between 
the peak frequency and the rms of all Lorentzians below ~20 Hz, 
i.e., when one Lorentzian shifted to a frequency at which an- 
other Lorentzian used to be, its rms was the same as the other 
Lorentzian had at this frequency. This relation could be inter- 
preted as a self-similarity, i.e., one Lorentzian profile replaces 
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Fig. 15. Dependence of the rms of the Lorentzian profiles on the photon index (left) and on the peak frequencies (right). The open 
squares correspond to the first Lorentzian and the filled diamonds to the second one. 



another one when it shifts to its frequency and originates from 
the same process or the same mode. Figure Q3] does not show 
this behavior. Owing to the different energy dependence of the 
Lorentzians it is obvious that this relation depends on the con- 
sidered energy channels. It is nevertheless possible that one 
Lorentzian replaces another one at its frequency, because when 
Li shifts to the position of L2 other properties of the source also 
change. For example V2 is about 5 Hz for T\ « 2.0, whereas 
vi reaches this frequency at T\ ~ 2.7. This change in the X-ray 
spectrum implies changes in the emitting region, possibly affect- 
ing the properties of the variability components. 



5.1.3. Time lags 

Previous time lag studies of CygX-1 showed typically short 
lags around 2-3 ms in the hard state, while enhanced lags 
up to 9 ms were observed during "failed state transitions". 
Pottschmidt et al. (2003) used the energy bands 2^-keV and 
8-13keV to obtain these values. When the source entered 
deeper into the soft stat e, the time lag returned to short values 
dPottschmidt et al.ll2000h . 

Our new data set with higher temporal resolution (RXTE 
orbit- wise extractions) allows us to extend these results (Fig.fTTl 
[T21) . In the hard state the time lag increases with softening X- 
ray spectrum consistent with previous observations, whereas af- 
ter the transition to the soft state the time lag behavior changes. 
Depending on the brightness, the time lag changes from values 
of up to 11 ms at low count rates to values of 2-3 ms at high 
count rates without appreciable softening of the source. Changes 
from long to short time lags occur from one RXTE orbit to the 
next. Figure QT| shows that the lag is short when the variability 
contribution of the power law is large, whereas the time lags of 
the Lorentzian variability are long depending on their peak fre- 
quency and thus on the spectral hardness. 

Time lags can originate from rep rocessing of radiation and 
propagation in the emitting region. iNowak et alJ dl999bl) ex- 



clude disk-intrinsic time delays as the sole reason for the lag, 
and discuss constr aints on properties o f the corona if the time 
lag originates in it. iKvlafis et al.l ([2008) present a jet model that 
can reprod uce the behavior of the time lags in the CygX-1 
hard state. Kordin g & Falckd (120041) show that time lags can 
be caused by a pivoting power law spectrum arising, e.g., from 
a jet/synchrotron model. In their work the sign of the time 
lag depends on the hardness/flux correlation ( Kordin g & Falckd 
2004, Fig. 9). Interestingly, the correlation visible in the lambda- 
shaped HIDs seems to change its sign coinciding with the change 
in time lag behavior, exactly at the transition: in the hard state 
the count rate increases with decreasing hardness, while it is the 
other way around in the SIMS. The transition is also where the 
behavior of the time lag changes. 

lUttlev et all d201ll) show that the time lags seen at soft 
X-rays (0.5-9 keV) are caused by the accretion disk. They 
point out that Comptonization models for lags have difficul- 
ties, whereas propagati on delays (see, e.g.. iKotov et al.l 12001: 
lArevalo & Uttlevll2006l) can explain their data, and that a sim- 
ilar mechanism can cause time lags at higher energies. 

The clear correlations between spectral and timing proper- 
ties put strong constraints on models explaining the emission of 
accreting black holes. 



5.2. Summary 

The main results of this paper can be summarized as follows: 

- CygX-1 was in an intermediate/soft state during 2005 
February 01-11. The photon index F\ reached values be- 
tween 2.0 and 2 .7 and was correla ted with F2 in the same 
way as found bv lWilmse t al. (2006). 

- The source was very variable during that period of time; 
the total PCA count rate varied between 1000 and 5000 



counts s per PCU. 
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- A full transition from the hard state to the soft state occurred 
in less than 2.5 hours. The two states can be clearly distin- 
guished in a hardness intensity diagram. 

- Spectral and timing parameters are strongly correlated. The 
peak frequencies of the Lorentzian profiles shift to higher 
values with increasing IV Those are identical in the low- 
and the high-energy band, whereas their rms depends on the 
energy. 

- In the HIMS the time lag increases with decreasing hardness. 
After the transition to the SIMS the time lag behaves differ- 
ently. At high count rates (>2500 counts s -1 PCLT 1 in this 
case) the time lag drops to about 2 ms and the contribution 
of the power law noise is large. At lower count rates the time 
lag is long, with values up to 1 1 ms. 
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